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In this paper, the size of bubbles formed through the breakup of a gaseous jet in a 
co-axial microfluidic device is derived. The gaseous jet surrounded by a co-flowing 
liquid stream breaks up into monodisperse microbubbles and the size of the bubbles 
is determined by the radius of the inner gas jet and the bubble formation frequency. 
We obtain the radius of the gas jet by solving the Navier-Stokes equations for low 
Reynolds number flows and by minimization of the dissipation energy. The prediction 
of the bubble size is based on the system's control parameters only, i.e. the inner gas 
flow rate Qi, the outer liquid flow rate Qo, and the tube radius R. For a very 
low gas-to-hquid flow rate ratio {Qi/Qo — > 0) the bubble radius scales as r^/R oc 
\/^Qi/Qo, independently of the inner to outer viscosity ratio r/i/r/o and of the type 
of the velocity profile in the gas, which can be either flat or parabolic, depending 
on whether high-molecular- weight surfactants cover the gas-liquid interface or not. 
However, in the case in which the gas velocity profiles are parabolic and the viscosity 
ratio is sufficiently low, i.e. rji/rjo <^ 1, the bubble diameter scales as rb oc {Qi/QoY i 
with ^ smaller than 1/2. 



PACS numbers: 47.55.db,47.61.Jd,47.15.Rq 
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FIG. 1. Schematic representation of the three main microfluidic geometries used for the formation 
of droplet and bubbles. The dispersed phase is injected (a) in a cross-flowing stream through a T- 
shaped junction, (b) in a co-flowing stream, and (c) in a focused stream imposed by the continuous 
phase. 

I. INTRODUCTION 

The controlled formation of monodisperse microbubbles and microdroplets at high pro- 
duction rates is important in many industrial and medical applications. For example, the 
food industry seeks new methods to generate en masse monodisperse droplets and bubbles 
to accurately control the density and structure of a wide variety of materials^ ^. In medical 
ultrasound imaging, microbubbles can be used as ultrasound contrast agents (UCA), where 
the resonance size of the bubble determines its acoustic respons^^. Commercially available 
UCA are produced using ultrasound-induced bubble formation methods, which results in 
a wide size distribution, i. e. UCA have a mean of 2 fim bubbles, but bubbles with a size 
between 1 /im and 20 fim exist in the population. Consequently, resonance occurs only for 
a small selection of bubbles. Improving the sensitivity in diagnostic imaging can thus be 
achieved by narrowing the size distribution so that more bubbles are at resonance size. 

The use of microfluidic technology for the generation of these accurately produced droplets 
and bubbles has received considerable attention recentlj^^. Various microfluidic geometries 
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are used, for example, T-shaped device^*^, co-flow, and flow-focusing geometrie^nHSl^ 
is shown schematically in fig. [TJ In the latter, the dispersed phase is focused by the outer 
continuous phase to enter a narrow channel (see fig. Droplets (or bubbles) are formed in 
the 'dripping' mode or 'jetting' mode depending on the system's parameters. These include 
the inner and outer volumetric flow rates, Qi and Qo, the channel dimensions, and the fluid 
parameters (viscosity and surface tension). When the outer continuous flow rate Qo is low, 
the droplets are generated in the dripping regime. In this regime, the dispersed phase enters 
the narrow channel and almost completely blocks the continuous phase. This leads to the 
formation of a neck, connecting the droplet to the inner feeding channel, that is gradually 
squeezed by the outer fluid until it breaks and a droplet is releasecP^. In Raven et al^ 
this mode of operation was used to generate a wide variety of foams. For increasing gas 
fraction (Qi/Qo) foams are produced that consist of separated bubbles, to bubbly flow, to 
dry bamboo foam. Droplet formation in the dripping regime is characterized by droplets 
with a size comparable to the size of the channel and its droplet production rate is typically 
low. 

When the outer liquid velocity is sufficiently high, such that the capillary number of 
the outer fluid Cao = rjou/'j > 1, with outer liquid viscosity rjo, velocity at the interface 
M, and surface tension 7, viscosity overcomes surface tension forces and an elongated jet is 
formecpS'^. This regime is referred to as the jetting regime. The cylindrical jet breaks up 
in equally-sized fragments driven by the classical mechanism of droplet formation through a 
Rayleigh-Plateau instability driven by surface tension forced^. High-throughput monodis- 
perse droplet formation in the jetting mode is of great value for industrial applications where 
a high-production rate is essential. 

Utada et al^ and Utada et al^ presented a simple description of the droplet size in 
a liquid-liquid co-flowing microfluidic system with a geometry similar to that considered 
here. In their work, the inner liquid viscosity rji is lower than the outer liquid viscosity rjo 
(rji/rjo ~ 0.1) and both the inner and outer flows are pressure driven; consequently, they 
exhibit Poiseuille-like velocity profiles (see fig. [2^). 

In this work we study microbubble formation through the breakup of a cylindrical jet in 
an axisymmetric gas-liquid co-flowing device formed by two coaxial tubes of circular section. 
Differently to the cases considered by Utada et all ^^ * ^^ l, we assume here an inner gas flow 
with negligible viscosity (r/i ^ rjo)- Two different situations are considered depending upon 
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FIG. 2. Velocity profiles for the inner and outer flow in a co- flowing stream, (a) Poiseuille velocity 
profile for both inner and outer flow for two co- flowing liquids, with rji < rjo. (b) The inner gas flow 
shows a flat velocity profile, whereas the outer flow shows a Poiseuille profile. The rigid interface 
allows for a discontinuity of shear stress at the gas-liquid surface. The red core represents the inner 
phase. 

whether the gas-liquid interface is free or rigidified by high-molecular-weight surfactants such 
as phospholipids, amphiphilic polymers, etc. These complex molecules are highly relevant 
in UCA microbubble formatiorP to avoid gas dissolution into the carrier liquid. Notice that 
the main difference between the two situations studied here is that high-molecular-weight 
surfactants naturally populate the gas-liquid interface and form a rigid interface that greatly 
affect the inner flow boundary conditions (see fig. |2|3). 

In this paper we obtain a complete description of the radius of the inner gaseous jet and 
bubble size solely based on: the gas and liquid flow rates, the size of the channel, and the 
liquid properties. We flrst describe droplet and bubble formation in the absence of high- 
molecular-weight surfactants, following the calculation already presented in Utada et al^. 
We do so in preparation for the key new flnding of this paper in which bubble formation 
with high-molecular-weight surfactants at the interface is described. We use minimization 
of energy dissipation to close the system of equations describing these speciflc condition^. 



II. DROPLET AND BUBBLE FORMATION FROM A LIQUID OR GAS 
JET 

In order to be able to highlight the analogies and differences between gas-liquid and 
liquid-liquid systems, we flrst repeat and elucidate the essence of the droplet-size calculation 
of Utada et al^ and then we extend this result to the case of bubbles. Let us consider 
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FIG. 3. Coordinate system for an axisymmetric unperturbed inner jet with radius f surrounded 
by a co-flowing liquid in a tube of radius R. Subscripts i and o refer to the inner and outer phase 
respectively. 

the flow of two immiscible fluids without surfactants in an axisymmetric co-flow device, 
as depicted in flg. These type of devices are typically fabricated by careful alignment 
of two coaxial capillaries: a tapered inner capillary for the supply of the dispersed phase 
and an outer capillary with radius R that delivers the continuous phase. When the outer 
liquid velocity is sufficiently high, such that viscous forces overcome surface tension forces 
(Cao > 1), an elongated liquid jet is formed. In the case of the dispersed phase is a liquid, 
the cylindrical jet ultimately develops undulations driven by surface tension that lead to 
the jet disruption into droplets with sizes comparable to the jet diameter. This breakup 
mechanism is known as Rayleigh breakupP^. The size of the droplets is = A^vrf^, where 
A* is the wavelength of the fastest growing disturbance and f the unperturbed jet radius. 

In this section we give an expression for the radius of the jet f as a function of the inner 
and outer flow rates, Qi and Qo respectively, the properties of the two fluids, and the outer 
capillary tube radius R. The coordinate system that represents the initial state of the jet is 
shown in flg. [3j 

The Reynolds number, as a measure of the relative importance of inertial forces to viscous 
forces, is typically low in microfluidics, hence the fluid flow can be described using the steady- 
state Stokes equations for low Reynolds number flows, 

= -Vp + riV^u, (1) 

with pressure gradient Vp, viscosity rj, and the velocity u. Note also that, with independence 
of the value of the Reynolds number, equation ([T| is still valid to describe strictly parallel 
streams, which is the case under consideration here. The velocity flelds are obtained under 
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the assumption of a no-slip boundary condition at the outer tube wall 



Uo{R) = 0, 



(2) 



continuity of velocity at the liquid-liquid interface 



Uo(r = Uiir 



(3) 



and continuity of shear stress 



dUn 
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(4) 



The unperturbed flow of the inner jet resembles a perfect cylinder with constant radius 
r = r. The capillary pressure, i.e. the pressure difference across the interface, is given by 
the Young-Laplace equation 

Pi-Po = i, (5) 



r 



with 7 the interfacial tension. Note that, since the radius of the cylinder is constant, the 
pressure gradient in both the inner and outer fluid are the same, and equal to Vp. 

Inserting the boundary conditions of eq. ([2]), ([s]), and ^ into eq. ([T]) and integrating 
gives an expression for the inner and outer liquid velocity profile 
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The velocity for the inner disperse phase and the outer continuous phase as a function of 
the radial coordinate are schematically shown in fig. [2^. 

The volumetric flow rate is found from the integration of the flux over the cross-sectional 
area of the jet as Q = 27r / ru{r)dr and, consequently, the flow rate ratio becomed^ 
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with X = r/R the dimensionless jet radius. This equation can of course be inverted, leading 
to 

, where X = — 
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FIG. 4. (a) Normalized radius of the inner jet as a function of the inner and outer volumetric 
flow rates for various viscosity ratios. The jet's radius decreases, while the maximum inner velocity 
increases, for decreasing inner-to-outer viscosity ratio. The radius of the jet scales as f oc (Qi/Qo)"- 
(b) Local slope a = dlog(f/i?)/dlog((5i/(5o)- In the limit of Qi/Qo — )• the size of the jet scales 
as f/R = 2^^/^ {Qi/Qo)^^"^ independent of the viscosity ratio. For decreasing viscosity ratio {i]i/r]o) 
a scahng a = 1/4 becomes pronounced. 

The dimensionless jet radius x as a function of the flow rate ratio (Qi/Qo) for various 
inner-to-outer viscosity ratios is depicted in fig. |4j It is shown here that, for Qi/Qo — ?■ 0, 
the size of the jet f/R = 2^^/^ {Qi/Qo)^^"^, independently of the viscosity ratio. The size 
of the jet is thus proportional to the square root of the flow rate ratio, a fact that was 
reported previously by various group^I^'^. Note, however, that for the relevant case of 
bubble formation, i.e. rji ^ rjo, the jet dimensionless radius is given by f/i? oc {Qi/Qo)" 
with Qi/Qo < 1 and l/4^a;^l/2fora large range of values of the flow rate ratio. Indeed, 
in fig.|4}D the crossover of the local slope a = dlog(f/_R)/dlog(Qi/Qo) from the 1/2-scaling to 
the 1/4-scaling is shown for various viscosity ratios. Only for the limiting case of an inviscid 



gas, say rjilrjo < 0.001, does the 1/4 scahng becomes pronounced. Note that a 1/4 scahng 
exponent is also observed in the co-flow device by de Castro et al."^^, but under the different 
conditions of a strongly pressure gradient in the entrance region of the outer capillary. 

For the case of drop formation, the cylindrical liquid jet is unstable against surface 
perturbations with a wavelength that exceeds the circumference of the jet (A > 27rf^)'2Sl. 
The fastest growing disturbance A* that leads to droplet pinch-off determines the droplet 's 
volume V = vrf^A*, with A* ~ 11. 2f (for rj^ = rjo^, hence the size of the droplets 
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Td ^ 2.03f ^ 1.44i? (Qi/Qo) • However, for the case of gas jets {r]Jr]o <^ 1), since the 
period of bubble formation is proportional to n R'^r{Qi/Qo), bubble volume is calculated 
through the mass balance 4/37rr^/6 oc tx R^Qi x f/Q^ ^ . This gives rj, oc {Qi/QoY, with 
the exponent /3 varying from the two limiting values /3 = 0.5 for Qj/Qo — ^ (a = 0.5) and 
/3 = 5/12 for rji/rjo — (a = 1/4), being these latter exponents flrstly reported by de Castro 
et al¥^. 

III. BUBBLE FORMATION FROM A 'HOLLOW JET' 

We now come to the speciflc case of UCA microbubble formation through the breakup 
of a 'hollow jet'. The co-flowing outer liquid contains high-molecular-weight molecules like 
surfactants, proteins, and amphiphilic polymers that self-organize at the gas-liquid interface 
to form a coating around the bubble. This coating, or shell, ensures a much longer bubble 
lifetime in contrast to uncoated bubbles, therefore improving their stability, an essential 
condition for medical imaging. 

In this section we give an expression for the size of the jet based on the inner gas and 
the outer liquid flow rates and under the assumption that the interface is rigidified by the 
surfactant solutiorP^. 

We solve the velocity flelds in the Stokes equations ([T| for an outer liquid flow rate with 
no-slip condition at the outer tube wall. We assume that the viscous stress imposed by the 
liquid is entirely balanced by the interface, hence that no stress is transmitted to the gas. 
This implies a homogeneous pressure distribution inside the jet, hence a flat velocity proflle 
(see flg. |2)d). 

The inner and outer velocity flelds respectively read 

Ui(r) = u = constant, (10) 
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FIG. 5. Energy dissipation e as a function of the radius of the inner gas jet n, for gas-to-hquid flow 
rate ratios Qi/Qo equal to 0.001, 0.01, 0.1, 1, and 10 (arrow indicating increasing gas fraction). 
For each given Qi and Qq there exists a minimum energy dissipation in the system, which defines 
the optimum radius of the inner jet x = f/R. 

with u the velocity at the interface, and 

Uo{r) = -— [R -r } + ^[u + — Vp . (11) 

4?7o \ogr/R \ Ar]o J 

Thus, the corresponding volumetric flow rates are given by 

Qi = vrf^M and (12) 

Qo RUl-x^)^ f 2 1 -2x2loga; 

— -Vp \ X — 1 + 



27r IGr^o V log a; 

Qi — 1 — 2x^ log X 



(13) 



vrx^ 4 log X 

with X = f/R. In this expression the pressure gradient (Vp) and the radius of the jet 
(x = f/R) are unknown. Experimentally, both are selected by the system once the gas and 
liquid flow rates are imposed. Here, it is possible to close the problem by minimizing the 
viscous energy dissipation 

/ M ^ ) dr. (14) 



2tx L V dr 



Inserting eqs. (11) and (13) into (14) and simplifying gives 



Qi J ?>{Qo/Qi-K{x)/2x^)^ 

(15) 



1 

x^ log X 



with K{x) = (x^ — 1 —2x^\ogx) /logx. Thus, the amount of energy that is dissipated 
into the system depends on both the inner and outer flow rate and the radius of the jet, 
as is demonstrated in fig. [5j The curves indicate the dissipation energy for fixed inner 
and outer volumetric flow rates and its minimum value marks the optimum radius x. The 
optimum radius as a function of the flow rates is found numerically by minimizing the energy 
dissipation function and is plotted in fig. [6] (solid hne). 

The gaseous jet is inherently unstabld^S and will break up in bubbles by means of either 
a capillary instability or as a consequence of a pressure drop in the gas. In the case of 
bubbles are formed as a consequence of a capillary instability, the breakup of the gas jet is 
driven by the fastest growing disturbance. A linear stability analysis predicts the optimum 
wavelength for the breakup of a gaseous jet surrounded by a liquid is slightly larger when 
compared to its inverted system — a liquid jet in air. The stability of a unconfined 'hollow 
jet' in ambient liquid was accomplished by Chandrasekhai'^Sl, who obtained that the size of 
the bubbles rb ~ 2.15f oc {Qi/QoY^"^ ■ In confined geometries as the ones considered here, 
the prefactor of this relation will be somewhat differentP^*^. In this case the jet breaks as 
a consequence of a pressure drop in the gas stream and the bubble formation frequency is 
proportional to T^rjQj' ^ . The bubble size can then be calculated from the mass balance 
rl oc R^fQi/Qo oc R^iQi/Qo)^^"^ and, thus, in this case, Vb/R oc {Qi/Qoy^"^- Therefore, 
independent of the source of the instability, the bubble size is proportional to the square 
root of the flow rate ratio. 

IV. DISCUSSION AND CONCLUSION 

In conclusion, we have expressed the radius of gas jets formed in axisymmetric co-flowing 
streams as a function the control parameters, i.e. the flow rates of both the continuous and 
the dispersed flows and the outer tube radius. The study has been divided in two parts, 
depending upon whether the gas-liquid interface contains high-molecular-weight surfactants 
to avoid the rapid dissolution of the gas into the carrier liquid. In the case where these 
complex molecules are not present it is found that, if the flow rate ratio is Qi/Qo — ?■ 0, the 
diameter of the bubbles is proportional to the square root of the flow rate ratio, independent 
of the viscosity ratio. However, for finite values of the flow rate ratio Qi/Qo < 1 and very 
low values of the viscosity ratio, i.e. rji/rjo <^ 1, the bubble radius scales as rb oc {Qi/QoY-, 
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FIG. 6. Optimum radius of the inner gas jet f as a function of the gas-to-hquid flow rate ratio 
Qi/Qo- The sohd hne represents the optimum radius that is found from minimization of energy 
dissipation and under the assumption of a flat velocity profile. Alternatively, the system is described 
by two pressure driven co- flowing fluids (ryi, 7/0)1 with parabolic velocity profiles for both the inner 
and outer phase for an air-water system (ryo = 10~^ Pa-s, rji = 1.8 x 10^^ Pa-s) (dashed line). For 
Q\/Qo — ^ 0, the slope is 1/2, indicating that the optimum radius scales with the square root of the 

1 /2 

flow rate ratio {f cc R [Qi/Qo) )• 

with (3 varying from the two limiting values (see fig. |4]) (3 = 0.5 for Qi/Qo and (3 = 5/12 
for rji/rjo — )■ 0. This latter exponent was reported by de Castro et al^. 

In the case that high-molecular-weight surfactants are added, we have found that the 
bubble radius is proportional to the square root of the flow rate ratio. Moreover, in flg. [6] 
a comparison is made between the gas jet radius predicted based on an inner parabolic 
velocity proflle (dashed line) and a flat velocity proflle (solid line) for an air-water system 
(with viscosity of air rji = 1.8 x 10~^ Pa-s and water rjo = 10^^ Pa-s). It is demonstrated that 
for small gas fractions Qi ^ Qo the influence of the inner flow conditions on the jet radius 
is marginal. 
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